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Abstract 

The effect of a strong magnetic field on the dominant neutrino emissivity in 
strange stars is investigated. In an ultra strong magnetic field, there exists an en- 
hanced neutrino emission because the charged particles are confined to the lowest 
Landau level. The results show that the neutrino emissivity is proportional to T'^ 
instead of conventional T", which implies more rapid cooling behaviors in mag- 
netars than usual stars 

PACS numbers: 26.60.+C, 97.60.Jd, I2.38.Mh, 97.60.Gb 

1 INTRODUCTION 

The theoretical possibility that strange quark matter, which made up of roughly equal 
numbers of up, down, and strange quarks, may be more stable than atomic nuclei 
(specifically iron, which is the most stable atomic nucleus) constitutes one of the most 
startling predictions of modern physics (DEI- If Witten's idea is correct, strange quark 
stars (strange stars) may exist (T'4l. Theoretical and observational researches on 
strange stars might provide a scientific basis to test what is the true ground state of 
the hadron. 

Initially a strange star cools mainly via the neutrino emission from its core. The 
(3 decay of quarks in quark matter is believed to be the important cooling mechanism 
in strange stars, which was usually investigated in absence of magnetic field in many 
worksri6"17V However, Observations of pulsars predict large surface magnetic field, 
typically of order lO^^G. Some magnetars are observed to have magnetic field of 
10^'* ^ lO^^G. Considering the flux throughout the stars as a simple trapped primor- 
dial flux, the internal magnetic field may go up to lO^^G or even more. In spite of the 
fact that we do not know yet any appropriate mechanism to produce more intense field, 
the scalar virial theorem (5) indeed allows the field magnitude to be as large as 10'^"G. 
In the case of super strong magnetic field such that' 2eB > pp{e)'^ J2Q'\, all electrons 
occupy the Landau ground state with electron spins pointing in the direction opposite 

*email address: zhxp@phy.ccnu.edu.cn 

'For electrons having a Fermi momentum of 50 Mev one obtains B > 2 X lO^'^G, for quarks such as d 
quark having a Fermi momentum of 300 Mev one obtains B <2x lO^^G. 
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to the magnetic field, and the charge neutrality now forces the degenerate quarks also to 
occupy the lowest Landau level even if the magnetic field is not large enough to break 
the Fermi surface. We also know that the quantization effect on the charged particles 
is extremely large when the magnetic field exceeds a critical field strength Be defined 
by equating the cyclotron energy ^ to to, (i denotes electrons and quark flavors) 1 191 . 
where the radius of cyclotron is —j^comvatibletodeBroaliewavelenatho fthevarticle f2^ . 

Be will respectively reach to ~ lO^'^, 10^^, 10^'* gauss for electrons, u(d)-quarks and 
s-quarks. Obviously, the effect on s-quark is small due to the large quark mass if the 
internal magnetic field is not over lO^^gauss. In our calculations, we consider the situ- 
ation that all charged particles involved in the processes confined in the lowest Landau 
level and it will give some instructions for general treatment regardless of the internal 
field in stars. In this situation we must use the exact solutions of Dirac equation. Emis- 
sion processes will be found to be evidently different from the field-free case. Many 
works recognized the magnetic effect on weak interactions in 90's of last century f6l[71 . 
In neutron stars the magnetic effect on the direct Urea, scattering and absorption pro- 
cesses was studied extensivelvllSl 1201 . We will show that the neutrino emissivity of 
ultra magnetized strange star is strongly dependent on the magnetic field. 

This paper is organized as follows. The exact solution of the Dirac equation for 
charged particles in the magnetic field is obtained in Sect. 2. The direct Urea processes 
dominating the neutrino emissivity in strong magnetized strange stars are calculated in 
Sect. 3. The related cooling curves are presented in Sect. 4. Finally we summarize our 
conclusion and discussion in Set. 5. 



2 EXACT SOLUTION OF DIRAC EQUATION FOR 
CHARGED PARTICLES IN MAGNETIC FIELD 

When an electromagnetic field with 4-potential is included, Dirac's equation is 
modified by the minimal coupling procedure of replacing by — qA^^, where q is 
the charge of the particle. The Dirac Hamiltonian is then replaced by 

if = a • (p - qA) + /3m + q4) (1) 

Here we set 4) = Q, just considering the charged particles in a constant, uniform mag- 
netic field in the z-direction. We choose the Landau gauge 

A = (0, Bx, 0) (2) 

and assume a trial wave function of the form 

'ijj{t, x) — f{x)exp{—ieet + iepyU + itpzz) (3) 

where e = ± is the sign of the energy whose magnitude is e. On inserting the trial 
solution (3) into the Dirac's equation in the form 

(iSt- ff)V'(i,x) = (4) 
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one can obtain the energy and the wave function. The energy is 

£ = + m2 + 2n\q\B 



(5) 



where n is the Landau level. It is degenerate and can be expressed as other quantum 
numbers such as the orbital part / and the spin quantum number s. 



n = l + -{l±s) 



(6) 



here the ± is on the contrary for the negative and positive charge respectively. One has 
q=\q\ior positive charged particles and their four solutions are 



v;(i,x) = 



for spin up, and 



exp[—ieet + iepyH + iep^z] 



( (e£ + m)In;p^ {x) \ 


\ -i^/2n\q\BIn-l■,p^{x) J 



(7) 



Vl(i,x) = 



exp[—ieet + iepyH + i^Pzz] 
y/2e{ee + m)LyL;^ 



( \ 

(eg + m)In -i;p, (x) 

iy/2n\q\BIn-p, (x) 

V -ePzIn-l;pAx) ) 



(8) 



for spin down cases. In Eqs. (7) and (8), Ly and are the length scales along the 
y and z directions. For negative particles q = — the wave functions are somewhat 
different as 



exp\—ie£t + iepyy + i(.pzz\ 
-y/2£g(e£, +m)LyLz 



exp[—ieet + iepyU + iep^z] 



( (ee + m)/„_i;p^(a;) \ 


epJn-i -p, {x) 
V i^2n\q\BIn-pAx) j 

(e £ + m)I n;p, (x) 
-i^2n\q\BIn-i-p, (x) 

-epzln;p, (x) J 



(9) 



(10) 



in above four Eqs. (7)-(10) 



^n;pz 



\q\b 



exp[-l\q\B{x - X ' 



^_^r.[VMBix-^)] (11) 



where Hn is the well-know Hermite polynomial. It is convenient to write 

4 = 



then Eq. (11) can be rewritten as 

, _f\<l\B 



(12) 



(13) 
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3 NEUTRINO EMISSIVITY 

The following weak processes are important in quark matter consisting of up, down 
and strange quarks, 

d— >u + e~+z>e, u + e~^d + Ue. (14) 

s^u + e~+Pe, u + e~^s + Uf.. (15) 

We will employ the Weinberg-Salam theory to calculate the neutrino emissivity 
of these reactions in strong magnetic field. The interaction Lagrangian density for the 
charged current reaction (14) and (15) maybe expressed as Ctnt = 
where G/ ~ 1.435 x 10^''^ erg cm~^ is the Fermi weak coupling constant and dc the 
Cabibbo angle. The lepton and nucleon charged weak currents are = ■i/'e7(ii(l ~ 
'y5)tp„ and = VuT^Cl - Ts)^^ respectively. 

when the u, d quarks and electrons are Landau quantized. The emissivity due to 
the antineutrino emission process in the presence of a uniform magnetic field B along 
z-axis is given by 

L,dvi, r L.dp, 



J-\q\BL^/2 J_ 



27r 



L^dpu. /-I^l^^^/^ ^^^^ 



J-oo 27r J_ 



ey 



\q\BL^/2 27r 



Aq\BL^/2 


Lyd/P^y 


J-\q\BL^/2 




r L.dpez 


r°° vd^pu 


J-oo 27r 


loo (27r)3 



xe^Wfif{ed){l - /(£e))(l - /(£„)) (16) 

here we consider all charged particles are in the lowest Landau level. The prefactor 
3 takes into account the d quark's flavour degenerate. The function /(£») denotes the 
Fermi-Dirac functions for the ith particle. The transition rate per unit volume due 
to the antineutrino process may be derived from Fermi's golden rule and is given by 
Wfi = (I Sfi \'^)/tV. Here t represents time and V = L^LyL^ the normahzation 
volume. The «S matrix element is given by 

Sfi = ^j d^X^.^^il - 75)V'i7e?„7^(l - 75)V'd (17) 

when all the charged particles in the lowest Landau level (n = 0), their wave function 
can be expressed very simply as 

eB 1 

V'e(x) = ( — y/*—^==exp[-ieSet + iepeyV + iepezz\exp{-^l)Ue,-i (18) 
cB 1 

V'd(x) = { — )'^^'^^==exp[-iesdt + iepdyV + iepdzz]exp{-£,j)Ud,-i (19) 

V ^y^z 

V'«(x) = (^^^1/4 _^^^ exp[-ieedt + iepdyV + iepdzz]exp{-^l)Uu,+i (20) 
V-i/.+ilx) = ij^j J exp[-ieeyt + iep^xX + iepuyV + iep^zz]U,y^+i (21) 

V J-ixJ-iyJ-iz 
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yLxLyL . 



(22) 



whereCe-(ei?)i/'(a;+|f 
The negative charged particles in the lowest Landau level spin down, while the positive 
charged particles spin up. So from Eqs. (7)-(10), we know that 



Ue.-l 



Ud. 



A/2e£e(e£e + ' 



\ 



u, + l 





e£e + me 




eed + md 


-f^Pdz I 





V 



Puz 





V2e£r(ee7 



/ e£y + 

ePi^z 



e£y + TO,, 

^Pux ^^Puy 
-<^Puz 



(23) 



(24) 



(25) 



(26) 



(27) 



Because only positive energy solution will contribute IIIOI (e = +1), using these wave 
functions it is straightforward to calculate the transition rate per unit volume which is 
given by (for convenient, we take p^x — Pv,y — Pv,z — 0.) 



6y/TTL,j;^Li 
5{pdy - Puy ~ Pey ~ PyVr)^{Pdz " Puz 



Pe 



-Pzv^) 



exp[ 



Puy 



^Ply 



Pey PuyPdy PuyPey PdyPey -i 



eB 



eB 



(28) 



eB 2eB 2eB eB 

here M is the invariant matrix element. Squaring it and summing over neutrino states, 
we get 



\M\' 



(Sd + md + Pdz)'^{£e + TOe - Pez)^(£u + + Puz)'^ 
4£d(£d + TOd)£„(£„ + mu)ee{ee + TOe) 

[(£i7, + m77, - Pzv^ + pIv^pIv} 

ei7,(£T7, + TO77J 



(29) 
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when the charged particle confined in the lowest Landau level, its energy is e = 
•\/p|~frn^. So if the particle's mass is zero, we have e — ±pz- Here we consider 
electron, u and d quarks' mass are zero and neglect the neutrino's momentum, there- 
fore only if 

Pdz>0, Pez,<Opuz>0 (30) 

the |Af p is not zero, it is 

|Mp = 16 (31) 

and note that the range of x can be fixed in a natural way 1 1 1 1 corresponding to 
the magnetic field = 1 /VeB. Considering all these factors, the neutrino emissivity 
can be written as 



'> ^6 r- / , "-Pdy / , "-P^v / , ' 

{2TT)°y/TT J_075^/2 J-y/2/3^/2 



Xfr,, „n -r, V-r-Tif ^"^ ^^''^ P^^vPdy PuyPey 

0[Pdy Pey Puy)exp[ cB 

PdyPey ^ / dpezdpdzdpuzd^pZzSu fd{l - fe){'^ ~ fu) 



SiSd - e„ - Ee - evJSiPdz - Puz -Pez ~ Pv^z) (32) 

Integrals over pdy, Puy and p^y can be performed numerically. The integrals over 
Pez, Pdz and puz are converted into integrals over dee, dsd and dsu respectively. The 
integral over neutrino solid angle is performed by using the zcomponent momentum 
conserving delta function. Finally using the standard procedure 1121 . the neutrino emis- 
sivity can be cast into 



TT^^^F^^ C 3 + 15C 5 dpey 



^■^(B) = , J.Z C(3) + 15C(5)]r^ / _ dp, 

Puy 5 (^Pey Puy ) Pcy Puy {,Pey ~t~ Puy ) 



dpuycxp 



^ eB 2eB 2eB eB 

(33) 



PuyPey . {.Pey ~l~ Puy^Pey 



eB eB 

Given the magnetic field, the integration in Eq. (33) can be calculated numerically and 
it varies about linearly according to the magnetic field. So the neutrino emissivity can 
be approximately written as 

e^{B) w 6.03 X IQ^^eBl^T^ (34) 

the neutrino emissivity of reaction (15) is similar to that if the magnetic field strength 
is higher than 10^^ gauss. 
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4 COOLING CURVES OF STRONG MAGNETIZED 
STRANGE STAR 



Because of the coulomb barrier, strange stars may have a similar thin crust to neutron 
starsO- We postulate that the crust of strange stars has the same structure with the 
same mass neutron stars. In our calculations, we take M = IAMq and R ~ lOfcm, so 
the corresponding central density of the the strange stars is pc — 6.65 x cm^^ . 

With regard to the relationship between the interior temperature T and the surface tem- 
perature Tg, we use the isothermal core approximation which assumes that the internal 
temperature T is constant within the stellar core with Pc > Pb — lO^^g CTO~^(Glen 
& Suthland 1983) and apply the result of Potekhin et al. (1997). We need discuss the 
heat capacity besides neutrino emissivity for the calculation of the cooling curves. We 
adopt the standard thermodynamic formula (Lifshitz and Pitaevskii, 1980). 

■^^ fd'p^is,-pM (35) 



where gj, Sj and pj stand for the particles' degeneracy, energy and momentum respec- 
tively, is the chemical potential, and fj is the Fermi-Dirac distribution. In field-free 
case, we have for an almost ideal, strongly degenerate ultra-relativistic gas 

Cj = — - — (36) 

In magnetic field, it reads 

<^. = f^/*«fe-M,)f (37) 

If we assume the charged particles to be in the lowest Landau level in the strong mag- 
netic field. We obtain 

gjeBT 

Cquark — ^ (38) 

6 

We can now make cooling simulation which is reduced to solve the equation of thermal 
balance 

where C, i-y and L^, denote respectively total heat capacity, neutrino and surface lumi- 
nosity, C = J2 J CjdV, L,y = J e^dV, = AnR^aTf with Stefan-Boltzmann const 
a and surface temperature Tg. Here we will compare two cases of the s-quark polarized 
and without polarization, which are showed in figure 1 and 2. In both the cases studied, 
the magnetic field results in more faster cooling compared to the the field-free case. 



C— = - (39) 



5 SUMMARY 

We have calculated the neutrino emissivity of strange stars with Ultra strong mag- 
netic field in some approximation and obtained an analytic formula. Our results show 
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that the neutrino emissivity is proportional to instead of as in the field-free 
case 11 61 . Bandyopadhyayl 18| had already calculated the neutrino emissivity of mag- 
netized strange stars, but they treated d and s quarks as free particles and they found 
that the relation between the neutrino emissivity and the temperature is the same as 
in the field-free case, with only a different coefficient. Further more, we find that the 
emissivity is strongly dependent on the magnetic field but independent of the electron 
fraction due to the extremely polarization. In the cases we studied, we show that the 
magnetic field accelerates the cooling of the strange stars. Finally, we emphasize that 
a general treatment must consider the contribution of different Landau levels of all in- 
volved charged particles because of the finite field strength. But we expect the neutrino 
emissivity is still proportional to T^, with only corrections to the coefficient, course, 
the analytical calculation is difficult and the numerical method should be applied. The 
future works will be necessary. 
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Figure 1: Cooling curves with the s-quark not polarized 
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Figure 2: Cooling curves with the s-quark polarized 
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